TECHNICAL MEMORANDUMS 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



No. 887 



COMPARISON OF THEORY WITH EXPERIMENT 
IN. THE PHENOMENON OE WING FLUTTER 
By P. Cicala 

L • Aerotecni ca , vol. XVIII, no. 4, April 1938 



Washington 
February. 1939 



J- 



f. J i- 



3 17 ^6 01308 7839 ! 



J 



NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



TECHNICAL MEMORANDUM NO. 887 



COMPARISON OF THEORY WITH EXPERIMENT 
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By P. Cicala 



The experimental data that are at present known with 
regard to the phenomenon of wing flutter permit essential- 
ly two classes of vibrations to be distinguished; namely, 
those for which it may be stated that the aerodynamic phe- 
nomenon follows the ordinary theory of potential motion, 
which theory assumes that the point of separation of the 
flow during the oscillation constantly coincides with the 
trailing edge of the wing; and vibrations arising from 
the irregularity in the aerodynamic phenomenon, whose pre- 
diction the present state of the theory does not permit. 
To the latter class belong those denoted by Studer as "de- 
tached" vibrations and perhaps also some types of slow os- 
cillation of the aileron.** 

A wing structure, in the normal range of angle of at- 
tack, may, up to a certain velocity, present oscillations 
of the second type. When the wing is brought to rather 
large positive or negative angles of attack, wind tunnel 
tests show a notable decrease in the critical velocity. 
In the wind tunnel the phenomenon presents itself as a 
serious one that would appear to justify the most pessi- 
mistic predictions. For a wing of strong camber free to 
undergo torsional displacements only and with sufficiently 
small friction at the supporting suspension, tests con- 
ducted at Turin showed rather low critical velocities in- 
dependent almost of the stiffness of the elastic suspen- 
sion system when negative angles of attack of only a few 



* "La teoria e I'esperienza nel fenomeno delle vibrazioni 
alari." L ' Aer o t ecni ca , vol. 18, no. 4, April 1938, pp. 
412-433. 

** There might also be distinguished forced vibrations 
(for example, oscillations of the plane of the tail due 
to the disturbed flow from the wing), but no sufficient 
data on these exist at present. 
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degrees were attained. The phenomenon, however, should 
be strongly affected "by the Reynolds Number so that from 
these tests there cannot actually be derived quantitative 
conclusions applicable to real structures. This aspect 
of the phenomenon i s at the present time being investi- 
gated at the laboratory at Zurich. 

The greater part of both the theoretical and experi- 
mental investigations has been directed to the study of 
oscillations of the first type. In the case of plane mo- 
tion (wing of infinite aspect ratio) the problem has been 
solved in all its aspects both as regards the aerodynamic 
phenomenon and as regards the conditions of stability of 
the system. To cite only a recent investigation, Kassnor 
(reference 1) has worked, out graphs, by the use of which, 
for any system of values of the fundamental- parameters, 
the solution of the problem may be obtained with suffix 
cient rapidity. When it is desired to consider the case 
ox the finite wing, tne study is considerably complicated, 
Important investigations have recently been conducted by 
the Aeronautical Ministry (reference 2). The equations 
which define the motion of the system are derived in a 
form analogous to that which is obtained by the use of the 
principle of virtual work (reference 3), the deformations 
being expressed by means of a linear combination of suit- 
able functions and the stability conditions discussed with 
the aid of Routh' s determinants. 

In all the studies so far developed, more or less 
arbitrary assumptions are made with regard to the aerody- 
namic actions. One of the most important of these assump- 
tions is that the forces which are developed at a section 
of the wing are independent of the state of motion of the 
adjacent sections. In order to remove such first approxi- 
mation assumptions, which seriously impair the precision 
of the methods of computation, there is required a knowl- 
edge of the aerodynamic forces and moments at the various 
sections of the wing corresponding to an assigned state of 
motion. . This problem has therefore been analyzed (refer- 
ence 4) by means of a proced-are which is founded on sim- 
plif iciations analogous to those which, in the study of 
the stationary motion of the finite wing (Prandtl vortex 
filament theory), lead to the familiar int egrodif f or ent i al 
equation of the circulation along the span. The present 
note gives the results of the numerical elaboration' of 
this theory. The coefficients of the aerodynamic forces 
and moments were determined for a wing of elliptic plan 
form for the two cases of aspect ratio 6 and 3, respec- 
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tively, and expressed in a . f orm .in. whi eh these coeffi- 
cients figure in the equation .of ..work. There, is also, in- 
dicated an . approximate procedure by the use of. which the 
resvilts may be extended to value.s of the aspect ratio and 
systems of deformation different from those considered in. 
the given computations. 

One result that appears from the- investigation is 
that the values of the coefficients of the aerodynamic ac- 
tions are only slightly influenced by the form of the, os- 
cillation. The computation carried out for one of these, 
coefficients (that Which gives the change in lift produced 
by the flectional (or flapping) displacement) gave almost- 
identical values for the cases of a deformation law repre- 
sented by a straight line and a second-degree parabola, 
respectively. This reult permits a considerable reduction 
in the labor of computation required for the applications. 
Another result of the analysis developed is that the coef- 
ficients, even for wings of normal aspect ratio, have val- 
ues deviating considerably from those given by the. theory 
of the two-dimensional or plane motion. This could have 
been predicted by the fact that in the vibrational phenom- 
enon the sections which undergo the greatest displacements 
are those at the wing tips most susceptible to the tip ef- 
fect. It is found therefore that the experimental investi- 
gations in which tnere tend to be realized the conditions 
of plane motion are valuable as confirming the theory, but 
do not furnish quantitatively valid data for the predic- 
tion of the critical velocities of actual wings which are 
of rather low aspect ratios. 

Also, in experimental investigations, the aerodynamic 
side of the problem has not generally been isolated and 
clarified. Normally the tests are carried out on a model 
wing in the wind tunnel by determining the. velocity at 
which the system begins to flutter and the corresponding 
frequency of the oscillation. Investigations of this kind 
have been conducted and others are being conducted in al- 
most all principal aerodynamic laboratories.. The diffi- 
culty which this method of investigation presents is that 
the results on the model are applicable to the wing if they 
present the same nondi men si onal ratios'. of mass arid stiff- 
ness and the same relative positions of' the. center of grav- 
ity and the elastic axis along the = chord. Thus, in order 
to realize in the wind tunnel test in a complete manner 
the conditions under which the effective structure is sub- 
jected, it would be necessary that 1 the model and the actual 
wing present at corresponding sections the same values of 
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the above nondimensio'nal parameters. On account of the 
extraordinary expansion in the extent of the investiga- 
tions required in this. case and the considerable difficul- 
ty encountered in the construction of models which are to 
have -the same- mass .relations as the actual wing, it seems 
more convenient to determine the aerodynamic actions "by 
diroct tests. The coefficients thus experimentally ob- 
tained substituted in the equations of stabi li ty ■ then per- 
mit the determination of the critical velocity and the 
other characteristics of the pehnomcnon for any system of 
values of the fundamental parameters, In addition to be- 
in,^ able thus to cover by a single series of tests the- 
range of -variation of the magnitudes which define the phe- 
nomenon, there is also eliminated the necessity of the 
construction of models which possess the assigned mass re- 
lations, Natural ly , in order that the results be applica- 
ble, it is necessary that the aerodynamic coefficients be 
independent of the Reynolds Number. If this is not the 
case, none of the critical velocity determinations thus 
far conducted on models would give results of practical 
value , 

On the basis of the above considerations there have 
been undertaken at the Aeronautics Laboratory at Turin 
direct measurements' of the aerodynamic actions on an os- 
cillating wing. The tests- so far conducted described in 
the notes of references 5 and 6 had as their essential 
object the examination of the operation of apparatus de- 
signed for this measurement and these tests will be re- 
peated and extended so as to cover a greater field of in- 
vestigation. The values experimentally obtained for the 
aerodynamic coefficients are in good agreement with the 
theory of oscillatory motion of the wing of finite span 
and show a clear deviation from the values obtained by 
the theory of plane motion. These de t er mi na.t i on s have 
permitted a oheck, given at the end of the present note, 
on the measurements carried out at the laboratory on the 
free oscillation of a model wing with variable stiffness. 
The agreement between the computation and experiments is 
satisfactory. However, various points must be clarified 
before being able to proceed with safety to the applica- 
tions. V7hen sufficient data has become available on the 
aerodynamic actions in the case of a wing with or without 
aileron and also with regard to the effects of the Reynolds 
Number, the mean angle of attack of the oscillation, and 
the compressibility of the medium; then the theory of wing 
vibration will be able to furnish reliable results for 
practical application. 
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Collection of Principal Symbols 
The principal symbols are the following: 
V, the wind velocity 
CI, angular frequency 
x, coordinate parallel to V 
L, local chord 
IT, mean chord (area/span) 

x cos 6, distance of a point of a profile section 

from the center point 

b , wing span 

z — - ^ cos t, coordinate normal to V 

A = — , aspect ratio 
L 

QL 

uu = — - reduced local frequency 

tt' = reduced mean frequency 

,u' = — 



+b /s 



- y -b/a 



Since the motion considered is harmonic, there is 
adopted the complex notation. All the magnitudes variable 
with time (displacements, vertical velocities, circula- 
tions, and lift and moment coefficients) are measured as 
deviations from a mean state, the exact knowledge of which 
is of no interest to the problem. These deviations are 
simply indicated by means of their complex amplitudes, sup- 
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pressing the factor e^*-. - •. The conditions relative to the 
mean state are studied separately by means of the equations 
of the stationary motion. The superposition of the effects 
is permitted because of the linearity of the fundamental 
e qua t i o n s . 

The flectional displacement is denoted by cp and the 
torsional displacement (rotation of the profile about the 
aerodynamic center) by . It is assumed that there ex- 
ists no phase difference between the displacements of the 
various sections. In every case it is possible to reduce 
the problem to these conditions by splitting up the dis- 
placements into two components in quadrature with each 
other and examining separately the two simple motions. 
The functions cp and \^ may therefore be considered as 
real. 

On account of the linearity of the fundamental rela- 
tions, the variation of lift produced by the flectional 
displacement is proportional to the displacement itself. 
Breaking up the constant of pr opo r t i onali ty i which in gen- 
eral is complex* into real and imaginary parts, we: write 
as for the plane motion (reference 5) 

c p = tt (a 1 *cp/L + ia^Gcp/v) (la) 

The nondimen si onal real coefficients a x * and a 2 * 

depend on the ratio QL/lT as for the plane motion, on 
the plan form of the wing and on the law of variation of 
the displacements. The latter vary along the span as a 
consequence of the fact that the actions on one section 
depend also on the state of motion of the contiguous sec- 
tions. 

Analogously, for the variations in lift produced by 
the torsional motion we write 

c p = rr(a 3 *\i/ + ia 4 * LQ\|//V) (lb) 

and for the moment about the aerodynamic center produced 
by the flectional and torsional displacements, respective- 
ly: 

Ca "= TTfb^cp/L + ib 2 *^cp/v) (1c) 



= n(b 3 *\|/ + ib 4 * Ln^/V) 
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.. . . The coeff icients which represent the actions:'' in phase 

(a j.*., a 3 * , , bx* , ./and. ..... b 3 *.).. ... contain th§,._lner t ia ..ef f ects 

of the air mass surrounding the. .wing. . This inertia action, 
Which remains independent of the velocity, .'is conveniently 
separated in the computations of the aerodynamic actions.^ 
If. the mass characteristics of. the wing are determined by 
a dynamic procedure, the actions are included in the meas- 
urement.... Making use of the results of tho theory of the 
wing of infinite 'aspect ratio (reference 5), we. separate 
the inertia terms (in u>2) and write 



»1* 


= a x «* 


+ 


<.u 2 




= a « * 

3 




a> 2 /4 


b * 
l 


= b x '* 


+ 


a) a /4 


b 3 * 


= b ' * 

3 




3 w 3 /32 



For the practical application, it is of greater interest 
to know the mean values of the aerodynamic coefficients 
along the span than the local coefficients and we define 
the former by the relations. 

■\ 

a i ,J b = / b a 1 »*qp a dz 

a 2 /- b cp 3 Ldz = a s * cp 2 Ldz 



a 3 1 •'b C9 ^ Ldz = •'V a 3 '*cpx|/Ldz 
a 4 / b cp\J/L 3 dz = a 4 *cp\|/ L 2 dz 
' J^cp \|/ Ldz = /-J, b x 1 *cp \J/Ldz 



> 



l3 2 / b cp\|/L a dz = /-jj' b 2 * cp \|/ L a d z 
b 3 *7 b \|/ a L 8 dB = / v b R «*\|/ a L 3 dz 



'b 3 



b 4 L 3 dz = /- D b 4 * \|/ I? dz 



(3) 
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It is easy to show that" the coefficients , a 2 , , . . 

t>4 defined by equation (3) satisfy the following condi- 
tion, namely, that the virtual work which the aerodynamic 
actions in phase and in quadrature with the flectional or 
torsional displacements perform against these displace- 
ments may be calculated by suppo sing • that the aerodynamic 
coefficients are constant for the entire span' and have' the 
values a 1 ',...b 4 given "by equation (3). A knowledge of 
these coefficients is therefore sufficient for the discus- 
sion of the stability of the wing structure by the method 
of virtual work (reference 3). It will be understood that 
the determination of these coefficients requires in each 
case the complete solution of the pro'blem for the assigned 
wing plan form and displacements. Since these character- 
istics have no great effect on the values of the mean coef- 
ficients, the solution of the problem for each particular 
case is not necessary. 



THE THEORY OF THE WING OE FINITE ASPECT RATIO 
IN NONSTATIONARY MOTION 



The problem of the oscillatory motion of a wing of ' 
finite span may be attacked without excessive complication 
when there are extended to it the simplifications assumed 
in the study of the stationary motion of the finite wing 
according to the vortex filament theory. 

Let c (fig. 1) be the wing contour in plan form. A 
vortex element a, which for the wing of infinite aspect 
ratio extends indefinitely in a direction normal to V, 
deviates in the case of the finite wing at a certain point 
and bends in the direction of the stream. In order to 
pass from the first to the second configuration we add to 
the vortex a the angle vortex h. If, in the computation 
of the induced velocity at a point 0, we substitute for 
the angle vortex h a parallel one having its vertex in 
Q, at the depth of 0, we arrive in the case of steady 
motion at the ordinary theory of the vortex filament. If 
the motion is not steady, there exist vortices with axis 
normal to V also in the wake behind the wing and there 
thus exist also angle vortices for the finite wing in the 
wake. The relations given in' this note are arrived at if 
the angle vortices whose vertices fall within the contour 
of the wing are displaced in the manner described for the 
steady motion and those of the wake are advanced by the 
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amount d, (fig... 2). , ..which represents the displacement 
ascribed to -these at the - instant at^which.. they •leave the 
wing. 

The sides of "the angle vortices perpendicular to V 
serve to eliminate the effect of the prolongation of the 
lifting vortices,, which are cut. off by : the section through 
the point 0 and which are considered as extended indef- 
initely on the two sides. The angle vortices having their 
vertices at the right of the section are therefore indef- 
initely extended . toward the right while those with their., 
vertices at the left extend indefinitely toward the left. 

As has teen said the effects of the. angle vortices 
are computed after displacements have been impressed upon 
them. This is equivalent in,each case to the addition of 
a transverse horseshoe vortex of which the inducing effect 
is neglected. ' (See reference 4.) The theory, is therefore 
approximate inasmuch as an entire system of vortices is 
neglected. In the case of steady motion, the theory leads 
to. r the integrodif f erential equation of Prandtl. Since the 
latter theory in ordinary cases leads to sufficiently sat- 
isfactory results, it may be expected that the extension 
to nonstationary motion leads to results of some impor-r 
tance, Iri order to arrive at more accurate results, it is 
possible, having solved the problem by means of the expres- 
sions given below, to determine the induced velocities of 
the neglected vortex system and treat these by the same 
procedure as the initial velocities and so on successively, 
Using the method of iteration.- 

In order to consider the effects of the nonstationary 
state of motion, we shall adopt, as for the plane problem, 
the method of Sirnbaum, which consists in breaking up the 
total circulation normal to the direction of V into a 
part that corresponds to the aerodynamic actions (bound 
circulation) and a part which- is shod from the wing in 
correspondence to the local variations in intensity and is 
carried along by the stream. (free circulation). 

lor the case of plane motion, it is known that the 
vertical (downwash) velocity distribution along the chord 
of the prof ile : represented by- one of the functions 



r 1 = v (-I + cos 

T n = V cos in £ (n = 2, 3. . , ) 
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corresponds to the distr ibuti on of the. hound circulation 

given respectively by the expressions 

r x = V ^2 sin £ - cot ^ - iiu sin i cos 6 - iu) sin £^ 

r = 7 { 2 sin n « — sin(n + l)d + sin(n - 1)*) 

ii v n + 1 n-1 / 



It is also known that the total circulation associated 
with the section corresponding to an arbitrary T n is al- 
ways zero, so that there are no free vortices in the wake. 
These asumptions lead to the conclusion that, for the finite 
wing, if we assume as valid the same simplifications that 
lie at the' basis of the vortex filament theory, the results 
holding for the plane motion may immediately be extended 
and it may be stated that,, if along the chords of the vari- 
ous wing sections, the vertical velocity distribution is 



expressed by means of 



7L 



or a linear combination .of them, 



the distribution of the' bound circulation will be repre- 
sented by the corresponding linear combination of r , as 

if. each section of the wing worked with an infinite aspect 
ratio without being influenced by the adjacent sections.. 



W 

"n 

of 



is not complete, however. It is 
the case of the rectilinear profile, 



except f z is zero 



bound circu- 



The series of 
sufficient to think 

for which every f n except f 1 is zero. The 

lation may therefore be expressed by means of T only if 
the vertical velocity correspond to the rear neutral point 
(cos -6 = -l/2) always remains zero, as results from the 
■expression W- .• 



In order to complete the series of 
the bound vortex distribution 

(2) 



r = 2 -a 

o 



H. 



+ i H, 



TIT 



n ' 



we consider 



co't — + iuusin 
2 



0 



in which A is a function of z and 
function of parameter uu . Denoting by 
culation associated with the section, we 



1 the Hankel 
K the total 
hav e 



c ir- 



K = tt A L/[X 
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where . ^ . , = f (h q ( 3 ) - i 

Since the total circulation associated with the various 
sections does not remain constant, it will be necessary to 
consider the wake vortices. With the simplifications pre- 
viously given, the s-cheme of figure 2 is arrived at in 
which is represented the configuration of the system of in- 
ducing vortices in the plane of the wing. In computing the 
velocity at a point 0, there are to "be considered: 

a.) The bound and free vortices a and I correspond- 
ing to the state of motion of the section through 0, which 
with the expression assumed for r x (see reference 7) in- 
duce along the chord a constant velocity given by 

-A = -u, K/rf L 

b) The angl.e vortices m whose vertices lie on the 
normal to V through 0 with density of circulation 
dK/dz, The velocity induced by all of these vortices ' is 



,■ ' d K d z 

I d z z - z i 

■ b 

c) The angle vortices h external to the wing. In 
the wake the vertices of these elements are distributed 
with density 

d z 

The velocity induced by one of these is 



/ 1 1 



4 TT \Z 



- + - / 2 + 

! Xj J (Z - Z x ) 



3 
1 / 



where x x and z - z x give the coordinates of the vertex 
of h with respect to 0. We introduce the function 



OO 



_ i i 1 + _1_ / 1 + 1 \ e ~ (XI 1 x x ^ 

./ \z-Zi x x V (z-z^ 2 x i 1 
o 
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The induced velocity of the vortex system c) is ex- 
pressed "by 



4 TT ,/ 



d K 
F d 

d z 



b 

Equating the induced velocity of the vortex system to 
W Q , which the stream should possess on account of its "be- 
ing tangent to the profile, there is obtained the integro- 
differential equation, which defines the function Z: 

u, 2 1 /' / 1 \ d K 

_ w = - - / ( tu'F)— -d z (4) 

o L 4 y v z - z x / d z 



If there is therefore given the vertical velocity dis- 
tribution expressed by means of the series 

W(x, z) = W Q (z) + W^z) (- + cos * ~\+ f 2 (z)cos 2 ■b + ... 

the total circulation K may be computed by (4) and the 
bound circulation distribution obtained by. the series 



2 (2) d 

r= f ,y rry- cot, - + i (ju sini V W^/V + W a r a /V +. 

tt L + i H 0 < 8 > 2 / (5) 

In the case of a nondef ormable section, the summation is 
reduced to the first two terms. 

From the bound circulation, there are calculated the 
aerodynamic actions bearing in mind that the difference in 
pressure between the two faces of the section is given by 
the relation 



p = P V r 

2 

Denoting the lift by c p L V dz and the moment 

of the aerodynamic actions about the aerodynamic center of 
an elementary strip by c p L 2 V 2 d z and making use of 

J? 

the expression (5) for T. there are found for the coeffi- 
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cients c and c , complex functions of z and of the 
state of motion of the system, the following expressions: 

c L V = / Td x = K + i tt (JUL(W 3 - W x ) /4 (S) 



P ./ 

-L 



1 + L / L 



c 

m F 



L 2 V = 7 r ( x + - )d x = i K L u) p./8 + 
•/_L 

tt W lL a (l - if-)/8 - tt W a L 2 (l - i^)/B - TT W 3 L 2 /32 

(a) 



where 

Hi 



Hj_ + i H Q d 



For practical application; it is necessary, first of 
all, to calculate the function F, which, with 

X = Qx/Y 

z = n ( z - z x ) / v 

may be written in the form 



(i) - /' + + 4 

, / V x z y x 2 z 2 ' 



1 X d X 



This* function is determined for possible values of 
the parameter Z. If Z is negative F is put in (4) 
with sign changed. This result immediately follows by 
considering the configurations of the vortices h to the 
right and left of the section through point 0 . Since 
this function is fundamental for practical applications, 
it was determined with particular care. For values of the 
parameter less than 4, there was made use of the approxi- 
mate expression 
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n Z e dX/ X 3 X 3 1.03 X 4 ' 

F(Z) = .-■■/ (i + ; > 

./ X + Z V 2(X + Z) 8(X + Z) 3 2(X + Z) 4 , 

' o 




For values above- 4', there was used the expression ob- 
tained by integration by parts 

On 

F(Z) = -i/Z + 1/2Z 3 - / 2. __±iE|2 - Z(3X 2 + 2Z S )(X 2 + Z 2 ) 

./ Q X3 L . J 

The integral in the second member was calculated nu- 
merically by extending the integration up to a limi;t such 
that tiie remainder, which could easily be estimated, re- 
mained less than the predetermined limits of approximation. 
The values obtained are given in table I. 

The function |j, of the parameter u; , also fundamen- 
tal for the computation, is given by the graphs of figure 
3 in which p,' and 1 are respectively the real and 

imaginary parts of the function. 



APPLICATION TO All ELLIPTIC WING 



The theory given above was applied to a wing with el- 
liptic plan form. In this case: 

L = 4 L ^i£_£/TT 

" ~Q - - 

We express the • total circulation by means of the sum- 
mation 

X = tt L V Z E n sin n i (7) 



where K n is the numerical complex constant. 
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1- 
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X 


? 4. 
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1 
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X 
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. o oy 
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JL 
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07? 
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4 
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3 0 


• \J DO — 
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1 


5 


. 778- 
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i 
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04ft 
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?74 


X 
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. 489- 
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i 


3 8 


o:^ o 
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x 
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.427- 
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i 
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1 
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.376- 


. 725 
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1 
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i 
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00c 
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X 
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i 
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i 




.022- 


.200 


i 


1.5 
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.567 


i 
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.020- 


.192 


i 


1.6 
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.540 


i 


5.4 


. 018- 


.185 


X 


1.7 


.176- 
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Equation (4) now "becomes 



-W Q sin t/V = 



TT |J, 



Z K n sin n £ + 



TT 



2nK n sinn I + 



4 



2 X 



+ i 



i - sin I / b F d £ S n K n cos n £ 



In the computation, only even terms of the series 
were retained so as to ohtain a velocity zero at the cen- 
ter section without having to add a supplementary equation, 
there being thus assumed .equal and opposite displacements 
for sections equidistant from the center section. More- 
over the state of motion of one half-wing has in our- case 
only a slight effect on the field of flow about the other 
half-wing. This is because in the ca.se of vibration the 
maximum displacements occur at the tip sections, which are 
therefore far removed or because the induced velocities of 
the angle vortices in the case of oscillatory motion de- 
crease with distance more rapidly than in the case of 
steady motion. 

Limiting the summation to the terms in K2 and K 4 , 
there were applied and. the valixes of ■ W 0 at the sections 
corresponding to z = C.7 and .z = b/2, and thus two 
linear equat i ons wer e obtained with complex coefficients in 
the constants Ka a-nd. K4 • 

It is convenient to consider the results separately 
for the various ■ component s . 

Lift produced by the flectional (f lapping) motion . - 
In this case there are constant velocities along each 
chord. Of the W n , only W Q . is different from zero' and 
given by 



where cp(z) is. the law which defines the f 1 ec t i onal . di s- 
placements along the span. 

Having determined the bound circulation in the manner 
explained above, it is necessary, for the computation of 
the mean aerodynamic coefficients, to know the quantity 



W 0 = i^cp 



Lrocp = / b c p L V 2 cp d z 



(8) 
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. , The real part of Leap represents the virtual work 
"by .... the aerodynami c actions in phase with the motion ; the 
imaginary part gives the w ork of the components in quad- 
rature. 

By (6) and (7), equation (8) "becomes 

L C pcp = PV/- b K|j, 1 cpdz =.tt P L V 3 / p, x cp d z S K n sin n t 

(9) 

3y (la), (2), and (3), equation (8) may also be written 

L cpcp = tt P V 3 (a 1 '/- b cp 2 d z +' / b cp 2 uu 3 d z + i^Sfifi / L cp 2 d z ) 

(10) 

Equating expressions (9) and (10) for !> corn , there is 
obtained an equation which gives the coefficients a 1 l and 
a 2 . 

The computation was carried oi.it assuming for the flec- 
tional displacements the law 

■P = - cos t I cos £ \ 

n 

In figure 4 is plotted the variation of the real part 
of the magnitude K/tt L V along the span for X = 6 and 
for the values of o> indicated on the figure. Figure 5 
gives the imaginary part. As may he seen from the curves, 
the total circulation varies appreciably with. the fre- 
quency of the oscillation. 

Figure 5 shows the distribution of lift produced by 
the flectional displacement considered. The curves give 
only the component in quadrature with the displacements 
(the ordinates of the cxirves give the values of the real 
part of Cp/rr). The component in phase has not been plot- 
ted, since it is of less importance. The reduced fre- 
quencies to which the curves refer are indicated on the 
curves. 

The values of the coefficient a 2 obtained for X = 
S and A = 3 are plotted in figure 7 as a function of (5. 
This coefficient decreases rather slowly with increase in 
the reduced' frequency, while in the case of plane motion 
it has a steep slope. The values obtained experimentally, 
indicated on the figure by small circles, have the same 
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general tendency as the theoretical values out lie some- 
what above the latter*. The measurements on this coeffi- 
cient will be repeated in order to ascertain whether the 
disagreement does not depend on the friction of the. oscil- 
lating system. 

The coefficient a. 1 t is rather small; the curve in 
figure 7 corresponds to A = 6. For practical applica- 
tions, this coefficient may be neglected and it may "be as- 
sumed, if the mass of the system is dynamically determined, 
that the lift is in quadrature with the displacement which 
generates it. Also in the test, the coefficient ^ 1 came 
out so small that it was- not measured. 

The computation was also performed assuming a linear 
law to represent the displacements; that is, putting 

cp = — cos i 
Q 

The values of the coefficients obtained in this case 
differ only slightly from those of the first case and were 
not plotted, since the curves would not be distinguishable. 

Moments about th e aerodynam ic center -produced by the 
flocti onal oscillat ion. - lor tiie- win,*, of infinite aspect 
ratio, the moment about the aerodynamic center produced by 
the flectional oscillation has no component in quadrature 
with the motion which generates it but is reduced to a 
pure inertia effect (b x = U) /4; b 2 = 0). The result may 
be held valid also for the finite wing. The component of 
the moment in quadrature with the flectional motion is 
equivalent to a displacement of the aerodynamic center by 
less than 5 percent of the chord 

/ K , \ 

( or i — - < 0.05 i . The component in phase except for the 
V i 3-2 : / 

inertia term is also negligible. 



*In regard to the comparison of the results of the tests 
conducted at Turin, it is necessary to observe that a cause 
of uncertainty exists in the fact that the aspect ratio of 
the model under the test conditions employed is not easy to 
define. The tests were- conducted in a free jet tunnel of 
rather small dimensions and the effect of tne jet boundary 
was rather difficult to estimate. The presence of a plane 
at the base of the model for the ptirpose of masking the un- 
derstructure complicates the phenomenon. 
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Lift produced by the -torsional oscillation. -- The .. 
focal axis is assumed rectilinear and there are deter- 
mined the aerodynaiai c • act i ons corresponding t o the tor- 
sional oscillations ahout this axis. In this case • \ 

-W = V ^ (1 + i.uu /2 - itu cos « ) 
and therefore 

-ir 0 = H(l + iui) : 
W x = i V \J/ (ju 
All the other W n are zero. 

Having determined the bound circulation, there was 
computed the work of the aerodynamic actions, due to the 
torsional motion, against the flectional displacements. 
By equations (lb), (2), and (3), this quantity is 

L ii-o = / t c p P L V 2 cp d z 

= tt P V 3 (a 3 ' /- 0 cp \|/ L d z - -i. / b u/ 3 .c? L d z + 

+ 2 — ; a t l d i) (11) 

■ v 

On the other hand, by (o) and (7), this quantity is equal 
to 

- 2 rr PV" S 

Hep = tt P L V cp (j^d z ZK n sin n £ + — - — / b iu 2 C9\|;L dz (12) 

Equating the .two, there are obtained a 3 ' and a 4 . 

In order to represent the angle of torsion, there is 
assumed the law 

v|/ = cos' t 

The values. of the coefficient a 3 1 for A . ;= 6 and 
X = 3 are shown on the curves of- figure 8. The coeffi- 
cient varies little with u5, which fact agrees with .the 
results o.f the test, indicated on the figure by small 
circles. 



20 



N.A.C.A. Technical' Memorandum No.. 38.7' 



Also the coefficient a 4 given"for X * 6 and \ = 
3 on figure 7 vari e s slightly, while, in the case of plane 
motion it decreases strongly toward small reduced frequen- 
cies. The experimental values indicated on the- figure "by- 
small circles are somewhat "below the theoretical. 

Mo ments about th e aero dynamic center due to the tor- 
si onal oscillation. - The work done by the aerodynamic ac- 
tions produced "by the torsional motion against the tor- 
sional displacements from (Id), (2 ) , and (3) is given "by 
the expression 

L ^ = rr p v 2 fb 3 '/ b * 2 L a d z - iL/ b u- 3 u/ 2 L 2 d z + ;*Vd z ) 

\ J 32 y / 

and from (6) and (7) by 

L\|n|/ = |tt P V 2 ( L / o L d z SK n sin n £ ■+ / b t .u (1 - i uj Mj/V'd z 

Equating these two expressions, there are computed the two 
coefficients b 3 ' and b 4 . The coefficient b 3 ', which 
for the infinite wing is zero, comes out as a small value 
for the finite wing. The coefficient b 4 , which for A = 
» is equal to 1.3, has for X = 5 and A = 3 the values 
given by the curves on figures 9 and 10. The experimental 
results are in good agreement with the theory. 

APPROXIMATE CALCULATION OP THE COEFFICIENTS 



Por systems of displacement different from those to 
which the previous curves refer, it is necessary to con- 
sider tne solution of equation (4). The computation, even 
if only approximate, is in every case very laborious: that 
of the velocities induced by an assigned system of vortices 
sufficiently simple for the case of steady motion must, 
for the nonstati onary ;aotion, be carried out by the method 
of graphical integration, and is complicated by the pres- 
ence of points of infinity in the function integrated. We 
shall therefore indicate an approximate procedure for the 
determination of the coefficients. The justification of 
the method lies essentially in the agreement of the re- 
sults with those obtained by using the more exact method 
indicated above. In the cases considered no practically 
appreciable difference was found between the results ob- 
tained by the two methods of computation. 
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From what has been said, whjat i„s "to be determined in 
each case is the work which. the aerodynamic forces due to 
an assigned system of displacements S perform when cer- 
tain .di sp lacements S' are impressed on the . system. The 
displacements under. consideration may be flectional, tor- 
sional, or those of the aileron, the condition required 
being that the motion of the- various sections be in phase 
If this condition is not satisfied, it is possible to re- 
duce to this condition by dividing the motion of the sys- 
tem into two components in quadrature. The displacements 
S and S' may be represented respectively by the real 
functions s(z) and s ! (z). With these functions known, 
there is first solved the equation 

s L = K_ - — [' (13) 



/ d. Zl 
b 



i " z 



The value of K 0 obtained by this equation represents 
the distribution of total • circulation , which under the 
conditions of steady motion would correspond to a dis- 
placement law given by s(z). The solution of the equa- 
tion may be found, for example, by following the procedure 
given in the paper by G-ebelein (Uber die Integralgleich- 
ung, etc., reference 8). 

If for each section there is computed the effect of 
the tip vortices, concentrating the circulation increase 
at the distance L/2 k on the two sides of the section, 
(13) becomes simply 

s L = K 0 ' (1 + k) (14) 



We define the numerical factor k (independent of 
z) by the condition that the virtual. work against the 
displacements s' of th« lifts corresponding to s, cal 
culated by the use of (14), have the exact value given by 
the solution of (13). That is, we write 

f h £ Q S' d-Z = £ 0 « B' d Z 



from which 



f . s s 1 L d z 

k(s, s') = -2 1 

A> K Q s' d z 
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Makinguse of this coefficient k also for the case 
of oscillatory motion, we concentrate the increase in total 
circulation at the same distance from' the section. as for 
the case of steady motion and equation (14) then "becomes 



-tt W 0 I = Z [(j, (u>) + k - i uu F(cu/k)] 



In order . to . compute the ..mean coefficients of the 
aerodynamic actions , we refer to a mean profile of chord 
IT. The total circulation will "be given by 



Z = 



-tt W 0 L 



[X (uO + k - i ai P (~/k) 
and therefore setting 



A = 



\x + k - i uu 



1 - A 1 - i A" = 



H 



(a) 



(2) (2) 



there will te obtained from, equations (6) 

i c V = -ff A (l - A' - i A" + U£) + iuJ (W a - 

TT V 2 / 



i c^v = iS^aw 0 + w,(i _ i_^\w 2 (l - iil)..- W 3 /4 



(15) 



The functions A' and A" are given in table I of 
reference 5. The values of \x are given in figure 3, and 
those of J in table I of the present note. The factor 
A, a function of u>, and through k, of the plan form 
of the wing and of the displacement laws, tends in the 
case of infinite aspect ratio (k = 0), to the value 
unity and in this case equations (15) furnish equations 
(1) of the above reference. 
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We observe that the value ' of k does not vary much 
on changing the system of displacements. For an elliptic 
wing of X = 6, there is obtained: 

For 8=8'= cos 1 k = 0.67 

and for s = s' = cost I c °s t ' k = 0.70 

L l 

The application of the method indicated to the case 
of flectional torsional oscillations is immediate. If only 
the first harmonic is considered, the laws of variation of 
the flecti onal and torsional displacements are not very dif- 
ferent. For the parameters k(cp, o) , k(\|/, c?) , k(cp, vj/) , and 
k(\|), \J/) , there will be sufficiently close values. Using 
a mean value of these, it will then be possible to determine 
the coefficients by means of the relations 

a + 2 i d) a 2 = -2 i 03 A (1 - X> _ i X" + i ffi/2) 

■ a 3 + 2 i o> a 4 = A(l + i £)) (1 - X' - i X" + i flu/2) 

b x + 2 i o> b 2 = A di 3 /4 f 

b 3 + 2 i dib 4 = i d>(l + i o>) A/8 + i o>(l + i (ju/4) /8 

which are readily deduced from equations (15). The exten- 
sion to systems with ailerons also is simple. 



CHECK OF THE FREE OSCILLATION TESTS 



From the measurements carried out, the mean values of 
the aerodynamic coefficients were derived and with these a 
check was obtained on the free oscillation tests conducted 
at the Aeronautical Laboratory at Turin (reference 9). 
For the computation, there has been assumed 

-a 3 = a 3 = 0. 55 



a 4 = 



0.38 
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- The- values of the other parameters of the structure 
are given in reference 3. The . cri t ical velocities are 
plotted in figure 11 as .a function of the ratio of flexural 
stiffness IU to torsional stiffness R^. . The experimen- 
tal values are indicated by the black circles. On the same 
figure are also plotted the values of the ratio of the fre- 
quency of oscillation corresponding to the critical veloc- 
ity to that of the pure torsion. The experimental values 
are in good agreement with the theoretical. 

»Ve denote by A the complex amplitude of the flee- 
tional displacement at the tip section of the model, by B 
the amplitude of the torsional oscillation and put A/BL = 
f + iq; ,f thus representing the component in phase and 
q the component in quadrature of the flectional motion 
with respect to the torsional. The curves obtained for 
these magnitudes are given in figure 12 as functions of the 
same parameter Hf/L 2 S t . If account is taken of the diffi- 
culty of such measurements, complicated by the fact that 
actually the condition of steady oscillations cannot be 
realized, it may be stated that the agreement of these val- 
ues with the theory may be considered as satisfactory. It 
is necessary also to observe that the two series of tests, 
namely, those of the free vibration and those of the aero- 
dynamic actions were carried out on entirely different pro- 
file sections. 

On the basis of .tne results obtained, it may be stat- 
ed that it is permissible, for an approximate calculation 
of the critical velocity, to apply the coefficients thus 
dot or mined . 



Translation by S. 
National Advisory 
for Aeronautics. 



Heiss, 
Commi ttee 
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Figs. 5,6,7 
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